Mean Field and Mean Ensemble Frequencies of a System of Coupled Oscillators 
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We investigate interacting phase oscillators whose mean field is at a different frequency from 
the mean or mode of their natural frequencies. The associated asymmetries lead to a macroscopic 
travelling wave. We show that the mean ensemble frequency of such systems differs from their 
entrainment frequency. In some scenarios these frequencies take values that, counter- intuitively, lie 
beyond the limits of the natural frequencies. The results indicate that a clear distinction should 
be drawn between the two variables describing the macroscopic dynamics of cooperative systems. 
This has important implications for real systems where a non-trivial distribution of parameters is 
common. 
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I. INTRODUCTION 

Systems consisting of large numbers of interacting 
units are common in science and nature, and have been 
the essential modeUing tools in physics, biology, chem- 
istry and social science [l[ . Here, the state of the whole 
system is characterized with macroscopic variables such 
as the temperature, magnetization and so on. The values 
of these depend on both the microscopic laws governing 
the dynamics of the units as well as their interaction. 
Generally, only macroscopic variables are accessible in 
experiments. Thus their precise definition and interpre- 
tation is needed. 

In the case of populations of weakly interacting oscil- 
lators, application of the phase approximation leads to 
the Kuramoto model (KM) for globally coupled phase 
oscillators 0. Although the model itself represents an 
idealized scenario, its analytical tractability makes it the 
prevailing approach in tackling a wide variety of irnpor- 
tant problems - from Josephson-j unction arrays [3| to 
brain dynamics under anaesthesia [3| and pedestrian in- 
duced oscillations in the Millenium bridge problem Q. 
This has led to many extensions of the basic model to al- 
low more realistic descriptions of actual systems, e.g. KM 
under influence of external flelds [sj , or with time- varying 
parameters |3] (for a review of generalizations and the 
problems they address see [8j and references therein). 

A fundamental feature of this model is that for a large 
enough coupling, synchronized behaviour emerges. De- 
pending on their inherent frequencies, some of the oscil- 
lators become locked, while the others continue to rotate 
asynchronously but with adjusted frequencies. The de- 
gree of the synchronization is usually characterized by 
some order parameter. For example, in the paradigmatic 
example of flashing fireflies [l| , this parameter would de- 
scribe the fraction of fireflies that flash in synchrony. 
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Since the building units of the KM are deflned by their 
natural frequencies, the macroscopic dynamics of the os- 
cillating system must be also characterized by some av- 
erage frequency. However, two quantities can be used: 
the effective frequency to which synchronized oscillators 
are locked and the average frequency of all the oscilla- 
tors, locked and unlocked, that belong to the observed 
system. The former represents the natural macroscopic 
frequency, whilst the latter is the microscopically aver- 
aged mean frequency. We shall call these respectively the 
mean field frequency and the mean ensemble frequency. 
Returning to the example of fireflies, the frequency of 
those flashing in synchrony - the mean fleld frequency - 
can generally deviate from the observed mean frequency 
of the whole population - the mean ensemble frequency. 
Similarly, some of the neurons in the brain are expected 
to be mutually locked to a certain frequency, whereas an 
electroencephalographic recording contains the mean of 
all neurons in some area, not just the synchronized. 

In general, there is no reason for these frequency def- 
initions to coincide. Still, not enough attention has 
been paid in formulating them for different parameters. 
Namely, due to the equality of the frequencies in the cases 
with symmetry that were mostly studied, they were used 
interchangeably and without veriflcation, even when they 
do differ (e.g. see [9]). However, here we consider sce- 
narios which most closely resemble the actual physical 
or natural phenomena; in particular, models with asym- 
metrically distributed frequencies, phase shifted coupling 
function, or asymmetric couplings of opposite sign. For 
them we show that these frequencies always differ and 
have non-trivial values. Hence, one should be extremely 
cautious when the measured frequency of a population 
is interpreted and then compared with the theoretical 
model. 

We begin with a formulation of the model and its group 
dynamics parameters. Section IIIII describes the station- 
ary solutions of the KM, whilst all possible scenarios with 
non-trivial mean field and mean ensemble frequencies are 
described in Section IIVI The summary of the work and 
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its implications are discussed in Section fVl 



II. FORMULATION 

The KM consists of phase oscillators running at arbi- 
trary intrinsic frequencies and coupled through the sine 
of their phase differences. In the case of heterogenous 
coupling strengths, the dynamics of the phase 9i of the 
ith oscillator has the form 
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9,), z = l,...,iV. (1) 



Here, Ki is the coupling strength of each oscillator and 
it is drawn from a probability distribution T{K). Simi- 
larly, the natural frequencies uJi are randomly distributed 
according to some g{uj). The tildes for the frequencies, 
phases and their distributions later are intentionally used 
for reasons to be explained below. Without loss of gen- 
erality, we assume g{ijj) to have a mean < Cji > centered 
at 0. Hereafter this frame of reference will be called nat- 
ural. Kuramoto introduced a complex order parameter 
for this model, defined as a centroid of the complex rep- 
resentation of the oscillators 
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It characterizes the macroscopic behaviour of the oscilla- 
tors, with r and "0 being the amplitude and the phase of 
the mean field respectively. The former shows the level 
of synchronization, while the latter gives the position of 
the peak in the distribution of phases. Generally, the 
long-term values for both can depend on time. Applying 
Eq. © to the governing equation ([T]), it is rewritten as 



9j = - K^r sin(6'j - 



(3) 



For infinitely large populations iV — > oo, a probability 
distribution function (PDF) f{9,uj,K,t) is defined, such 
that /+J fi9, Co, K, t) S = g{Q)T{K). Thus, the complex 
mean field Eq. ([2]) becomes 
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e"' f{9,Cj,K,t) d9 d£o dK. (4) 



For convenience, infinite limits in all further definite in- 
tegrals will be omitted. 

As a consequence of the conservation of the number of 
oscillators the evolution of the density function is gov- 
erned by a continuity equation 

Here the right hand side of Eq. ([3]) is used and the sine 
function is expressed with complex exponents. Moreover, 



since / is 27r periodic in 9 it allows a Fourier expansion 
and can be written as 
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where c.c. are the complex conjugates and = 

In the limit t oo, the ensemble described by ([5]) 
might settle into a stationary state for some rotating 
frame. We define: i) systems that have stationary so- 
lutions in some frame of reference i.e. the complex mean 
field and the distribution of phases rotate uniformly and 
they have a constant mean field after the initial transi- 
tions; and ii) systems that experience complex non-stable 
behavior, i.e. a time-varying amplitude of the order pa- 
rameter r{t). These definitions might differ from usual 
descriptions found elsewhere which regard as stationary 
only those solutions that are fixed in the natural refer- 
ence frame. In this work, our attention is focused on the 
ensembles with stationary solutions as described by i). 

For the case of identically coupled oscillators with uni- 
modal and symmetric distributions of their natural fre- 
quencies, above the critical coupling a phase locking of 
the oscillators takes place around the peak of g{u]) where 
the density of the oscillators is highest As a conse- 
quence of the symmetry, the group dynamics is station- 
ary and both mean frequencies are equal to the mode of 
g{u}), which in this case is also its mean value. 

Nevertheless, introducing multimodal or asymmetric 
5(0)), or distributed K leads to much richer dynam- 
ics. Thus, for multimodal g{6j) bistabilities and stand- 
ing waves emerge [lol - [l^ . A standing wave is a macro- 
scopic solution where neither f{9,ui,K,t), nor z{t) are 
stationary in any rotating frame. This further implies 
non-stationarity of r{t). Standing waves are also ob- 
served in systems with symmetrical bimodal distribution 
of natural frequencies, with the exact result for bifurca- 
tions between different states given in [T3| . 

Travelling waves (TW), as another peculiar group be- 
haviour, have also been observed for different parameter 
ranges within the same systems [13, El- In our analysis, 
a TW state is considered to be any solution character- 
ized by long-term stationarity of the mean field ampli- 
tude, whereas the frequency of the locking D, differs from 
the mean of the natural frequencies. In other words, the 
locking of synchronized oscillators is in a frame differ- 
ent from the natural. This also represents a stationary 
solution accordingto its definition above. 

A recent study shows the occurrence of TW in mod- 
els with positive and negative coupling strengths, and 
identifies so called conformists and contrarians. Simi- 
larly, a synchronization around a frequency that is dif- 
ferent from the mean or the peak of the distribution 
was reported for ensembles that have an asymmetric uni- 
modal distribution of natural frequencies IJ] . It is also 
worth mentioning here the Kuramoto-Sakaguchi model 
, where the phase shift is introduced into the coupling 
function, to allow synchronization at a frequency differ- 
ent from the mean of the natural frequencies. Hence, it 
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always leads to TW states. Additionally, the whole class 
of models of non-isochronous oscillators with constant 
shear can be reduced to this model [3, ■ 

Stationary solutions for fully symmetric populations 
have macroscopic frequencies that are equal to < >. 
In asymmetric scenarios on the other hand, the syn- 
chronized cluster experience non-trivial phase velocity 
[1, H; 13 • Thereafter, the focus of this work is inter- 
acting phase oscillators whose coherent behavior is char- 
acterized by a TW state, as defined earlier. Having cer- 
tain asymmetries, either in the frequencies, the coupling 
parameters, or in the coupling function itself, is a nec- 
essary condition for occurrence of this state. As a con- 
sequence, the influence of the unsynchronized oscillators 
to the entrainment frequency does not vanish 0, [l5| . 
Additionally, we show that these oscillators also cause the 
mean ensemble frequency to have a non-trivial value. It 
generally differs from the entrainment frequency, but also 
from the mean and the mode of the natural frequencies. 

Before proceeding with the analysis of all cases with 
TW state, let us define the macroscopic frequencies that 
we have already discussed. The mean field frequency rep- 
resents the velocity of the mean phase ip and is obtained 
from the time derivative of the complex mean field ^ 
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Taking into account that r and are both real, the same 
is true for their time-derivatives, so from Eqs. 1^ and (0 
one finally obtains following evolutions of the amplitude 
and the phase of the complex mean field 
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^^^ = ^T.^^s- - Vi)] cos(0~, - ^). (8) 



The expression for ip represents the velocity of the mean 
phase, i.e. the frequency of the synchronized oscillators. 
For the other frequency parameter - the mean frequency 
of the ensemble - used for characterizing these systems, 
its definition leads to 

1 ^ • 1 ^ 
/e„, ^-J2e,^- - K,r sinie, - V^)]. (9) 
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In the infinite limit, © is introduced into (j?]). By taking 
the time derivative and applying the substitution ([5]) , the 
evolution of the complex order parameter is obtained 

1 = 11 [zCof* + f (5 - ~z*f;)]mnK)d£:;dK. (10) 

Similarly, the mean frequency of the ensemble becomes 

fens= f I I ~e ]{e,u,K,t)dedCjdK. (ii) 



III. STATIONARY SOLUTIONS FOR THE 
PHASE DISTRIBUTION 

For a large class of problems, the long term coher- 
ent dynamics of an ensemble of phase oscillators is time- 
independent, as a consequence of the stationary distribu- 
tion of the phases. This implies existence of the station- 
ary solution of the continuity equation ([5]), which does 
not need to be in the natural reference frame. A recent 
work discusses generalized empirical stability condi- 
tions for these systems. 

From now on we consider that the ensemble has non- 
zero mean field, i.e. it is out of the incoherent state 
(a fully incoherent solution is also stationary), and that 
stationary solutions for the phase distribution exist in 
some planes of reference. This is true for the simplest 
possible scenario - unimodal symmetric frequency distri- 
bution and constant coupling strengths. The TW state is 
another possible scenario with this property, despite the 
fact that in this situation the velocity of the mean phase 
(or the mean field frequency as defined here) differs from 
the mean of the natural frequencies. 

In our analysis the mean field frequency is allowed to 
be nonzero despite assuming that < a) >= 0. Still, we 
work in the frame where f{9, uj, K, t) is stationary - the 
reference frame rotating with the frequency VL. Here the 
phases of the oscillators are 9 = — fit and the phase 
corresponding to the complex order parameter is ^ = ^p — 
fit. The distribution of the natural frequencies becomes 
g{uj) — g{uj + fl) with mean —Q. In the same frame "0 = 
can be assumed after an appropriate phase shift. 

For any ensemble, if the stationary solution of ([5]) ex- 
ists, then it exhibits two types of long-term behavior, 
depending on the size of |a) — f7| = |a;| relative to \Kr\ 
and to the sign of K. The oscillators with \uj\ < \Kr\ 
approach a stable fixed point defined implicitly by 



TT + arcsin ■ 



ii K >0 
ii K <0 



(12) 



These oscillators are called locked because they maintain 
constant phase difference, while rotating at frequency 
in the original frame. It is also assumed that synchro- 
nized oscillators with positive couplings have phases in 
the interval (— 7r/2, 7r/2), while those with negative are 
in (7r/2, 37r/2), since these are necessities for stable so- 
lutions of Eq. ([T^ Q. In contrast, the oscillators with 
\uji\ > \Kir\ rotate in a non-uniform manner. As ex- 
pected, the locked osciUators correspond to the center of 
g{uj) and the drifting oscillators correspond to the tails. 

For the synchronized oscillators the stationary distri- 
bution of the phases becomes 



g{Lo)T{K)S[e - encsh 
g{u:)T{K)S[0 -aicsin{j^)^7r 



ii K >0 
UK <0. 



Oscillators with frequencies in the interval |a;| > r\K\ are 
out of synchrony with the mean phase. Their stationary 
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distribution is obtained from the continuity equation ([5]) 
and the normahzation condition of f{6, uj, K, t), such that 



fUO,^^K,t)^,i^nK)^^^^^^. (14) 

Hence, the real and imaginary parts of the complex mean 
field definition, Eq. (|4]), become 



= 



cos 0f{e, UJ, K, t) de dio dK, (15) 
j J J sm0f{9,uj,K,t) d9 duj dK. (16) 



The latter is identified as the phase balance equation. 

The distribution of the phases for synchronized oscil- 
lators, Eq. (|13p . for each uj implies tt difference between 
the phases of the synchronized clusters with couplings of 
opposite sign. Nevertheless, this holds for the distribu- 
tion fs{uj, K,0,t) only if g{uj) is symmetric. This is not 
the case for the TW state, where even if the oscillators 
were symmetrically distributed in the natural frame, the 
symmetry would be broken when moving to the frame 
rotating with fi. Thereafter, the centroids of the phases 
will be shifted from the tt mutual distance, a phenomenon 
that was reported in [9,] . 

In the rotating reference frame the locked oscillators 
are frozen, i.e. by definition 6 = ijj ~ 0. Thus only the 
drifting ones need to be considered. This can be applied 
in the expression for the mean frequency of the ensemble, 
Eq. (fTTj) . which becomes (see Appendix |A|) 



fens^ [g{oj) - g{-oj)]r{K)^oj'^ - [KrYdujdK. 

J J\Kr\ 

(17) 

The definition ()17p is not restricted to any distributions 
of the natural frequencies and couplings. 

Going back to the rotating frame where the original 
g{Cj) has zero mean, frequency parameters of the TW 
state are characterized by £7 and fens^ where 
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As a consequence of ft can equal only if g{uj) 
is symmetric in the boundaries of the integral in p7)) 
[l3 |. This means that for any zero centered distribution 
of the natural frequencies g{uj) which is even in the 
intervals w > \rK\ (e.g. Lorentzian, Gaussian, etc) the 
function inside the integral is also even if 17 — 0. Thus 
the integrals cancel each other and fens = fens — ^ — 0- 



IV. THE TRAVELLING WAVE STATES 

In the following we discuss the possible scenarios that 
lead to the TW state for which we also obtain mean fre- 
quency parameters described earlier. 

All further analyses are carried out in the rotating 
frame of the entrainment frequency f2, such that ip = 
and natural frequencies are distributed according to 
g{Lu + n), where g{oj) has zero mean. However, for better 
clarity of the figures, frequencies in the examples with 
asymmetric g(tl') are depicted as seen from the original 
distribution described in the captions. In other words, 
compared to the results from the analysis, the plots are 
shifted by the means of the given distributions g{u}). 



A. Traveling waves in the Kuramoto model with 
contrarians and conformists 

First we focus on the TW solution described in Q, 
which actually inspired this work. The model shows re- 
semblance to sociophysical models of opinion formation 
and is also a continuation of the KM with distributed 
positive couplings [l9| . 

The distribution of the couplings is 

r{K) = {l-p)5{K-Ki) +p5{K-K2), 

where Ki < and K2 > 0, and p denotes the probability 
that a randomly chosen oscillator is a conformist, while 
q — 1 — p is the probability that a random oscillator 
is a contrarian. The natural frequencies follow a zero 
centered Lorentzian distribution with half-width 7 
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As stated in [9| , if the absolute coupling strength is higher 
for conformists than for the contrarians, then for some 
region in the parameters space 7 — p, the synchronized 
oscillators will experience a TW. This means that both 
peaks in the phase distribution uniformly rotate in same 
direction. The waves appear in symmetric pairs, with 
frequencies ±$7, since they result from the asymmetry 
in the coupling strengths. On the contrary, if the TW is 
due to the asymmetry in the natural frequencies or in the 
coupling function, the waves are not paired, as discussed 
later. 

Following the definitions of T{K) and g{oj), and using 
the substitution oj^ — {Ki/2r)'^ = W121 the integrals in 
Eq. ([T7|) are analytically solved, yielding 
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1-p 

^(72 + + i^2^2)2 _ 4/^2^2^2 _ f]2 + J^2^2 + ^2 



^(7^ + 1^2 + 1^2^2)2 _ 4K'ir^n^ + K^r^ + 72 ' 



(18) 
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FIG. 1. (color online) The amplitude r and frequency Q of 
the mean field, and the mean ensemble frequency /ens versus 
the ratio p. Theoretical results are given with a solid line 
(black) for r, a dotted line (red) for $1, and a dashed line 
(blue) for fens ■ Results from numerical simulations are shown 
with squares (black) for r, triangles (red) for Q., and diamonds 
(blue) for fens- Parameters: 7 = 0.05, Ki = —1 and K2 = 2. 



This expression can be straightforwardly generalized for 
multimodal-5 distributed coupling strengths. 

It is obvious that 17 = will imply /ens = 0. Hence 
only in the presence of a TW may the mean frequency of 
the ensemble differ from the mean phase velocity for this 
model. Similarly, for the TW state fens is non-zero and 
has opposite sign from il. Additionally it can be shown 
that the expression in the curly brackets is smaller than 
1/7, so that in the natural reference frame, \fens\ <\^\ 
always holds. This is also evident from the numerical 
results plotted in Fig. [TJ 

Let us now derive the expression for obtaining the fre- 
quency 17 of the TW, as seen in the natural reference 
frame. Ott and Antonsen, in their seminal work [2l| . 
showed that macroscopic evolution of large systems of 
coupled oscillators can be described by an explicit defi- 
nite set of nonlinear differential equations. They intro- 
duced an ansatz for the complex Fourier coefficients in 
Eq. © 



Uco.K,t)^[a{Co,K,t)r, 



(19) 



which exactly solves the governing equation ([5]), as long 
as a{Cj,K^t) evolves following the nonlinear equation 



da K 2 ~*\ 

-]- luja -\- — ( za — z ) — V 



(20) 



dt ' '"'^ ' 2 

When this ansatz is implemented in Eq.Q, the order 
parameter reduces to 



a*{uj,K,t)g{uj)r{K)dojdK. (21) 



I 

Note that applying this ansatz in pUj) . leads to a sim- 
plified expression for the evolution of the complex order 
parameter 

v&* + y (z - ra*^)]g{Cj)T{K)d(:}dK, (22) 

from where the frequency of entrainnient can be ob- 
tained. Similarly, substituting Eqs. ((S)) and ([6]) and the 
ansatz into Eq. (ITT|) transform the mean ensemble 
frequency to 



fens — 



[u;--(za- 



S*a*)]g{oj)T{K)d(:jdK. (23) 



We start by considering the low-dimensional evolution 
((20)) in the reference frame of the TW. Similar analysis, 
but in the natural frame was performed in Q- First 
the bimodal-(S and Lorentzian distributions for T{K) 
and g{uj) respectively, are substituted into the integral 
Eq. ([2T|) . such that it yields 



z* = {l-p)ai{n~i-f,Ki)+pa2{n-i-f,K2). (24) 
This is then used in Eq. ([20l) . which rewritten for both. 



q;i(17 — i"f, Ki) and a2(17 
dai,2 



17, K2)^ results in 



^ +(jl7 + 7)ai.2 + ^(^<2-^*) = 0. (25) 

where we omit the dependencies of a{uj, K^t). As pre- 
viously stated, we are interested in a stationary solu- 
tion of the TW state in the t ~^ co limit. This im- 
plies time-independent distribution of the phases in this 
limit, or from the ansatz (jl9p . time-independent ai,2 with 
dai,2/dt = 0. Further, similar to l9ij, complex order pa- 
rameters for each of the subpopulations, and the differ- 
ence between their phases are defined 
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"2, 



S ~ ipi — 1^2 ^ const. (26) 

In this way it is ensured that both synchronized pop- 
ulations, in-phase and antiphase, rotate with the same 
velocity 17 in the natural frame, and preserve constant 
phase difference. In the rotating frame of the TW tp = 0, 
so ■0 = can be set without loss of generality. Thus, from 
Eqs. ([25l - [26]) we obtain the following evolutions which 
describe a fixed point in the {ri,r2,ipi,tp2} space. 





-7^1 - -yii^i - 1)(P^2 


cos 5 + qri)] 


= 0, 


(27) 


r2 = 


-7''2 - -^[(''2 - 1)(P^2 


+ qri cos (5)] 


-0, 


(28) 


^1 = 


17 — -r^pr2 sin(5(r? + 1 
2ri 


= 0, 




(29) 


-02 = 


17 H -qri sinSirn + I 

2r2 


= 0. 




(30) 
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The low-dimensional parameters including n can be now 
obtained self-consistently. The steady states (|29j|30p re- 
sult from the constant angle difference between the peaks 
in the phase distribution, i.e. 

S ^-sm5[^pr2{rj + 1) + ^qnirj + 1)] - 0. 



Hence, when /to[z] 
use the expression 



2r2 

is applied to Eq. ([M)) . one can 



qri sin ijji = —pri sin ■02 



(31) 



and the definition (^5]) of 6 to obtain the values of ipi and 
■02, such that the system will be fully described. 

The equations (j27ll30p and all further numerical inte- 
grations are performed using a Runge-Kutta 4th order 
algorithm. Once we have the low-dimensional parame- 
ters, Eq. ([TOl) is applied to find the mean frequency of the 
ensemble. Finally, results are compared with the values 
for the order parameter r, the mean phase velocity of the 
ensemble and the mean frequency fens, obtained from 
the numerical simulations of the ensemble Eq. ([T]), using 
Eq. ([2]) and Eq. As in all later simulated scenarios, 
the number of oscillators was set to iV = 100000, the 
time step of the integration was 0.01. The simulations 
were running for 10^ time steps, with the initial 90% of 
each run discarded as possibly transient, while the rest 
were time averaged. The proportion of the conformists 
p is changed from to 1 at 100 equally spaced points, 
and the obtained results given in Fig. [T] fully confirm the 
theoretical analysis. 



B. Asymmetric unimodal frequency distribution 

Another case of the KM characterized with a station- 
ary solution, where the mean frequency of locked oscilla- 
tors differs from the mean of the oscillators' natural fre- 
quencies, are ensembles with asymmetric unimodal dis- 
tribution of the natural frequencies, and equal couplings. 
The asymmetric scenario is also more natural than the 
symmetric, because any imperfection in the system, how- 
ever small, can destroy the ideal symmetry. Nevertheless, 
mostly due to the analytical difficulties, this case has ob- 
tained little attention. It was first examined by Sakaguchi 
and Kuramoto [l^ and the self-consistent condition for 
the mean field frequency was obtained in 14] . 

It is known from 1J| , that the phase balance equation 
P6p for gild) asymmetric in the interval |a;| > rK , implies 
that the oscillators always lock to a frequency that differs 
from the mean of the natural frequencies. For the mean 
ensemble frequency Eq. (jl7p was numerically integrated 
for ensembles with triangular and lognormal frequency 
distributions and constant coupling strength. We found 
that it always gives fens = in the TW rotating frame. 
In other words, the mean ensemble frequency equals the 
mean of the natural frequencies, although it differs from 
the mean field frequency. The results seemed unexpected 
at first sight. However after careful analysis we realized 
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FIG. 2. (color online) The amplitude r and frequency Q, of 
the mean field, and the mean ensemble frequency /e„3 for 
unimodal asymmetric natural frequencies. Results from the 
numerical simulations are shown with squares (black) for r, 
triangles (red) for Q and diamonds (blue) for /ens. The lines 
(blue) are from theoretical results for fens and they also match 
the means of the frequencies' distributions - crosses (blue), 
(a-b) Triangular g{ij) in boundaries a = 1 and b — 10. (a) 
Coupling K = 4.2 and mode c G [0, 10]; (b) K £ [4, 10] and 
c — 9; (c) log-normal distribution of natural frequencies, with 
M = 0, cr G [0.25, 2] and coupling K = 4.2. The dashed line 
shows the modes of the distributions. 



that the derivations for the equation of balance Eq. ([2 
as given in (l^ , in the frame rotating with £7 indeed leads 
to the Eq. ([TT]). Namely, applying the derivations in the 
Appendix El Eq. (HH) becomes 



= 



\Kr\ 



giLu).yuj^ - {Krfduj 



\Kr\ 



g{0j)^0J^ - [KvYduj = n + fens. (32) 



Thus, one can conclude that for the KM with asymmetric 
frequency distribution, the mean ensemble frequency is 
always equal to the mean of the natural frequencies. 

Next, a triangular distribution with limits a and 5, and 
peak at c was explored. In the scenario shown in Fig. [5] 
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(a) the peak is distributed in the interval [a,b], while in 
Fig. [5] (b) the coupling strength is increasing from 4 to 
10 for fixed c. Similarly, for log-normal frequencies 

fi is fixed to 0, while a is logarithmically distributed, and 
the results are given in Fig. [2] (b). The mean ensem- 
ble frequency values also match theoretical values for the 
means of the given distributions, {a+b+c) /3 and e'^+'^ 
respectively. 

For the triangular distribution shown in Fig. [2] (a) the 
mode is c, while for (b) is fixed. So for the log-normal 
g{ui) given in Fig. [2] (b) the mode is 1 for /i = and 

2 

cr — ?> and then exponentially decreases to with e'*"'^ 
as shown. Hence, the presented results show that the 
mean field frequency is always between the mode and 
the mean of the distribution of natural frequencies, and 
reaches the second only when all oscillators become syn- 
chronized. That is to say, by increasing the coupling 
strength for a given frequency distribution, the propor- 
tion of synchronized oscillators is also increased. Accord- 
ingly the value of r2 moves closer to the mean of the dis- 
tribution, until it eventually reaches it for r — >■ 1. For 
unbounded g{uj) the last can only occur when K ^ oo, 
while in the case of bounded natural frequencies, for some 
value of rK all oscillators will be entrained. This can 
also be deduced from Eq. ([5^ . Namely, higher synchro- 
nization implies a smaller region for the integral on the 
right hand side. At the same time we assume g{u}) to be 
decreasing left and right from the mode, meaning that 
higher r leads to smaller value of fens which eventually 
becomes 0, implying O = fens- Furthermore, the value 
under the square root in the same integral also decreases 
with increasing r. These are confirmed in Fig. [2] (b-c), 
where we see that for smaller couplings and hence for 
smaller mean field amplitudes, is closer to the peak of 
the distribution and approaches fens for larger coupling. 

C. Asymmetric multimodal frequency distribution 

The KM with multimodal asymmetric distribution of 
natural frequencies is another candidate for the mean 
field behavior described by a TW. Nevertheless, due to 
difhculties that arise in the mathematical analysis of this 
model, it was never fully solved, nor has a thorough dy- 
namical analysis of possible macroscopic solutions been 
performed. Still, following the analysis of the symmetric 
scenario and qualitative descriptions of the dynamics 
in the asymmetric case given in j2| and [l^ , some conclu- 
sions can be drawn. Namely, Kuramoto, in his seminal 
work [2]; discusses how transition from incoherence to 
mutual synchronization might be modified when the os- 
cillators' natural frequencies are bimodally distributed. 
For sufhciently large coupling strength, he assumed that 
the clusters of synchronized oscillators "will eventually be 
entrained to each other to form a single giant oscillator" . 



For smaller coupling strengths compared to the distance 
between peaks, he envisaged that the synchronized nu- 
clei would be at the peaks of g{uj). Although some of 
the transitions between different states described in Q 
for the symmetric case were shown to be wrong , the 
description of the mean field dynamics of a partially syn- 
chronized state is indeed correct. 

Hence, for a symmetrical ^(w) the ensemble should be 
either partially synchronized for large enough coupling 
compared to the peaks' distance, or synchronized clusters 
should exist near the peaks 0, [l^ . For asymmetri- 
cally and bimodally distributed frequencies the partial 
synchronization will be characterized by TW, whereas 
standing waves and other complex collective rhythms ap- 
pear for smaller couplings. 

Thereafter, for the partially synchronized state in the 
asymmetric bimodal case, one might expect a TW state 
to occur. As a consequence of having a single synchro- 
nized cluster, the dynamic of the system is similar to the 
models with unimodal asymmetric distribution described 
in Section IIVBI As a result the same conclusion drawn 
about the asymmetric unimodal scenario, also holds for 
the TW solution of the asymmetric bimodal case. This 
means that the mean frequency of the ensemble is equal 
to the mean of g(a)), and as such differs from the fre- 
quency of the TW experienced by the mean field. Hence, 
in the TW frame, fens = This can be seen in Fig. [31 
where results from numerical simulations confirmed the 
theoretically expected values for fens and Q. 

As an example we analyze the dynamics of a popula- 
tion with bimodal Lorentzian distribution of frequencies 

~(~\ (1 --P)7i Pl2 

The subdistributions are peaked at ^1 and /i2 with the 
half-widths 71 and 72 respectively, where p is the pro- 
portion of the oscillators belonging to the second sub- 
distribution. The full low-dimensional dynamics for this 
system can be obtained using the ansatz Eq. (jl9[) in a 
similar manner as in Section |IV3 The integral (gl]) for 
the given T{K) and g{Lu) yields 

z* = (1 -p)aiipi - -ii,K,t) +pa2ini - ji,K,t). 

Hence, substituting ai 2 in Eq. (j20p . the low-dimensional 
dynamics are described by two ODEs 

K 

ai,2 = -(« Ml, 2 + 71, 2)"!, 2 - y (2* + zaia)- (33) 

In all the plots in Fig. |3] the values of the coupling 
strength are large enough to induce TW instead of stand- 
ing wave states, i.e. one instead of two synchronized 
clusters. For the Fig. [3] (a) both peaks are of equal 
distance from the zero. This leads to zero mean for any 
width of the subdistributions and following previous dis- 
cussion and Eq. (j32p it implies fens — 0- For similar devi- 
ations of subdistributions the entrainment frequency will 
be near zero, but closer to the narrower peak, since more 
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oscillators will asymptotically reach the peak of the sec- 
ond subdistribution in the limit case 71 — >■ oo. If in the 
same limit case p — > also holds, then with the same 
reasoning very few entrained oscillators will be close to 
the second peak. At the same time the mean frequency 
of the distribution and of the ensemble will be near the 
first peak. Similarly, for p = 0.25 and p = 0.75 as in 
Fig. [2](b-c), the mean of the frequencies is on 1/4 or 3/4 
of the distance between the peaks, respectively, while fl 
reaches the narrower mode. 

The generality of the described analysis for any bi- 
modal distribution is shown in Fig. [3] (d), where the 
frequency distributions are Gaussian. Although the low- 
dimensional dynamics in a form similar to Eq. p3p can- 
not be obtained for this case, the results from the numer- 
ical simulations of the ensemble, Eq. ([1]), are in line with 
the previous discussion concerning Lorentizan g{(jj). 
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FIG. 3. (color online) The amplitude r and frequency fl of the 
mean field, and the mean ensemble frequency fens versus the 
frequency width of the first subpopulation 71. (a-d) Results 
from the simulations, Eq. ([T]), are shown with squares (black) 
for r, triangles (red) for O and diamonds (blue) for /e ns . { a-c ) 
The low-dimensional dynamics, Eq. (|33|l , are shown with solid 
lines (black) for r, dotted lines (red) for Q and dashed lines 
(blue) for fens- (a-c) Bimodal Lorentzian g{ijj). Parameters: 
71 € [0,10] 72 = 0.75, m = -1, H2 = 1, K = 5, and (a) 
p — 0.5, (b) p = 0.25 and (c) p — 0.75. (d) Bimodal Gaussian 
g{u}) with 71 G [0, 10] 72 = 0.75, fii = —1.5, fi2 = 1, coupling 
K = 4.25 and p = 0.5. 



oscillators from that subdistribution will be entrained. 
However, for a large deviation of the first subdistribu- 
tion compared to the second one, a small number of the 
oscillators belonging to it will be entrained. These oscil- 
lators will be almost equally distributed on both sides of 
the second peak. Thus, the frequency of the entrained 



D. Contrarians and conformists with asymmetric 
unimodal frequency distribution 

The analysis naturally continues with the KM in the 
presence of both previously analyzed conditions that are 
required for TW - contrarians and conformists, with 
asymmetric distribution of natural frequencies. The 
phase locking in this scenario happens in frames other 
than the natural for all coherent solutions. But as it will 
be shown, due to the distributed couplings, the mean en- 
semble frequency also has non-trivial values - it is not 
always equal to < J) >. As another consequence of the 
asymmetry, the sign of the TW can no longer be expected 
to be random, as is the case for the contrarians and con- 
formists over symmetrically distributed frequencies Q. 
However, the analysis of this case is largely more com- 
plicated compared to previous ones and we have taken 
some examples only to show the main points. 

The phase balance equation for distributed K in 
this case does not lead to a simple expression as was 
Eq. ([221) for equal couplings. Still, following the same 
procedure, Eqs. (fT3l[T4)) are substituted into Eq. ([TG]) and 
applying the derivations from Appendix the balance 
of phases becomes 







TiK) 
rK 



dK{ / ujg{uj)duj 



a;[.g(w) - g{-uj)yuj2 - {Kr)^dLj} 

\Kr\ 

r{K) 



rK 



-dK [-^ + I{K)]. 



(34) 



This expression self-consistently gives the frequency of 
the entrainment, whilst it no longer implies equality of 
fens and < a) >, as for constant coupling parameters. 
Results from numerical simulations shown in Fig. |4] con- 
firm the non-trivial nature of both global frequency pa- 
rameters. The most interesting phenomena is that VL 
and fens for some parameters can have values that not 



only differ from either the mode and the mean of the 
natural frequencies, but that are also outside of the in- 
terval between them as was case for equal K. Moreover, 
they can even be outside of the boundaries of their dis- 
tribution. For example for the triangular g{Cj) g [0, 1] 
depicted in Fig. |4] (a), the frequency of the entrainment 
can be up to ~ —1 or to ^ 2, and similar happens with 
fens ■ These values are clearly out of the region of support 
of <D, which is shaded in the same plot. Note that these 
on first sight counterintuitive results are a consequence 
of the very high coupling strengths with opposite si_gn, 
and they are still within the general boundaries of fens 
and SI. Namely, from Eqs. (fill) and Q, it is clear that 

d) — max |iir| < /ens < cj 4- max |if I, (35) 

while from the fact that the entrainment frequency can- 
not be out of the limits [min(6'), max(6')] it follows 

min(a)) — max \K\ < < max(cj) + max \K\. (36) 

Of course, these are broad limits and only for bounded 
g{Cj) do the boundaries for f2 not reach ±cx). 

Hence, for bounded distributions, as the triangular, 
and for large enough couplings of opposite signs, the only 
way the balance equation (|34)) can be satisfied, is with the 
emergence of a TW with a large enough value, so that the 
integral I{K) in the same equation will be non- vanishing. 
Contrary, if I{K) — for the shown example, i.e. 

Co-n<\rK\, V w for which g(c;;) > 0, (37) 

then the phase balance becomes 

Q = {\-p)/Ki+p/K2. (38) 

If ([38| holds, then can have any value that still obeys 
condition (|37p . so each simulation would give a different 
TW within these limits. Another consequence from con- 
dition ([57]) when applied to Eq. (IT7)) is that fens ~ 0. 

If ((38|) is not satisfied, then the oscillators rearrange, 
rendering the condition (|37|) also invalid. The result is 
the appearance of non-zero I{K) in Eq. ([M]) to impose 
the balance of the phases. There may be more than one 
value for O which impose this balance and if that was the 
case one might expect multiple solutions or even hystere- 
sis behavior. In this work however, we only numerically 
confirm that for bounded natural frequencies, as in the 
example in Fig. |4] (a), SI, when observed in the natural 
frame, can have two stable values with opposite signs. In 
other words, numerical realizations could lead to any of 
the two stable values. This is illustrated in Fig. |4] (a) 
where one such realization is depicted, whilst in different 
realizations the values on both sides of the mean and the 
mode oi g{ijj) appeared for p around (0.35,0.7). 

For the same triangular natural frequencies, if p > 
0.65, then a state similar to the tt state in the symmetric 
model in Section IIV Al is observed. Here, the equality 
of SI and fens follows from the bounded distribution of 
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FIG. 4. (color online) The amplitude r and frequency fl of the 
mean field, and the mean ensemble frequency fens for asym- 
metric ensembles with distributed coupling strengths. Nu- 
merically obtained results, Eq. ([TJ, are shown with squares 
(black) for r, triangles (red) for and diamonds (blue) for 
fens ■ (a-c) The dashed lines (blue) are theoretically predicted 
results for /ens, Eq. (|17p . Horizontal dashed and dotted lines 
are the modes and the means of g{oj) respectively, and its do- 
main is shaded for (a), (a) Triangular 5(a)) within a = and 
6=1 and peak at c = 0.8; bimodal-J T{K) with Ki = —4, 
K2 = 8 and p £ [0.34, 1]. (b-c) Log-normal g{Cj) with jj, — 
and bimodal-5 r{K). (b) a = 1.1, Ki = -10, K2 = 20; and 
(c) a = 1, p = 0.55, Ki e [-6, 13], K2 = 6. 

the oscillators' natural frequencies, which all become en- 
trained. Similarly the integral I{K) vanishes, and the 
requirement for the phase balance holds only if SI = 0. 
In this way only the first integral of the equation p4p sur- 
vives and it gives < uj >, which we set to be (although 
in Fig. |4] the plots are for non-zero mean frequency dis- 
tributions, such that the plot (a) has mean 0.8). 

If the natural frequencies are unbounded then I{K) ^ 
0, meaning that a TW always emerges. For the simula- 
tions we have performed, as seen in Figs. |4] (b-c), the sys- 
tem was always setting the same value for n that solves 
Eq. (j34l) for the given parameter range. Thus, for the log- 
normally distributed w, SI is always positive, although the 



10 



existence of another stable TW that also fulfills the phase 
balance is not excluded. The entrainment of all oscilla- 
tors is never achieved in this scenario, and the mean fre- 
quencies will never reach < cD >. Still they will approach 
this value when the number of entrained oscillators is in- 
creased, either by increasing the number of conformists. 
Fig. U (b) , or by changing the coupling strength of one of 
the groups, as shown in Fig. U (c). The same plots also 
show that even when the couplings are all positive, the 
group frequencies behave qualitatively differently from 
the case with constant K. This characteristic becomes 
less obvious once the modes are closer to each other. To 
check the generality of the analysis, we also explored the 
case when K is bimodal Gaussian distributed. The re- 
sults did not show any qualitative difference from the 
case with bimodal-(5. 

It is expected that for multimodal couplings, cluster- 
ing emerges and plays a crucial role in defining the values 
of fe.nft and Q. A similar phenomenon was also reported 
in [l9| for distributed positive and unimodal K, but as 
a finite-size effect only. That is to say, the distribution 
of the locked oscillators Eq. ([T^ for multimodal-^ ^{K) 
will have peaks on different angles ipi — arcsin for 
each mode Ki . Further evidence for this explanation can 
be seen in numerical simulations performed over different 
T{K) and g{uj). Fig. [5] presents obtained phase distribu- 
tions. As can be seen three qualitatively different types 
were identified. 

Firstly for very large coupling strengths compared to 
the width of g{u}), and for bimodal T{K) with modes 
having different signs, there are two separated peaks 
corresponding to the contrarians and conformists. The 
distance between them is smaller than tt and the ob- 
served TW is similar as in Section IIV A[ although the 
asymmetry in g{uj) causes asymmetry in the distribution 
f{d,u!,K,t). This is shown in Fig. [5] (a), and also in ex- 
amples (a-b) of Fig. m in the range for p around 0.4—0.65. 

The second scenario differs from the first only in the 
sense that the absolute values of the coupling strengths 
are not much bigger than the width of the natural fre- 
quencies. As expected, this corresponds to the tt state 
from but the form of g{uj) will be mapped onto the 
peaks which are separated by tt. The asymmetry addi- 
tionally will influence the values of the fens and fl. How- 
ever, the latter will always be in the interval between 
the mode and the mean of g{uj). This means that qual- 
itatively this corresponds to the case with constant K 
described in Section IIV Bl but with non-zero values for 
fens because of the phase balance equation ([M]) . 

Lastly, if the couplings distribution have positive peaks 
only, depending on their distance apart, the phases will 
have either one or two close peaks, while still keeping the 
form of g{uj). Hence, the case with one peak will be the 
same as having unimodally distributed K and numeri- 
cal results suggest that fens ~ (w) as for constant K, 
although this is not immediately obvious from Eqs. (jl7p 
and Finally, if f{9,Ld,K,t) is no longer unimodal. 




FIG. 5. (color online) PDF of phases for different states of 
the ensemble with log-normally distributed natural frequen- 
cies and bimodal-(5 distributed coupling strengths. Oscillators 
with coupling Ki are shown with a dashed line (red), those 
with K2 are following the dotted line (blue,) while the joint 
distribution is shown with a solid line (black). Parameters: 
/i = 0, a = 1 and p = 0.5. (a) Ki ^ - W K2 = 20. Two peaks 
at a distance smaller than n and a TW occurs, (b) Ki — —4 
K2 ~ 8. Two peaks at a distance tt. (c) Ki = 2 K2 = 3. The 
peaks merge into a single peak in the joint distribution, (d) 
Ki = 2 K2 = 6. Two closely separated peaks. 

then fens has values that differ from < uj >. 

The number and the nature of the emerging clusters 
is not so straightforward for non-d F{K) distributions. 
Still, from the PDF of the phases it is clear that for neg- 
ative modes in r{K), the peaks in the PDF corresponding 
to those oscillators do not have the form of g{uj), but it 
is more symmetric. For positive modes, the peaks keep 
mapping the distribution of the natural frequencies. This 
is also an interesting peculiarity that requires further at- 
tention. 

The bottom line is that different and on first sight 
counterintuitive values for the entrainment and mean en- 
semble frequency are obtained even for all positive but 
multimodally distributed couplings, as long as the natu- 
ral frequencies are nonsymmetric. Nevertheless, a thor- 
ough analysis is still needed for this general case of the 
KM, with the emphasis put on the entrainment at differ- 
ent clusters and their influence to the group behavior. 

E. Phase shifted coupling function 

The flnal model leading to TW of the macroscopic pa- 
rameters being analyzed here is actually the earliest gen- 
eralization of the KM [1^ . This model was introduced to 
allow for entrainment in frames other than the natural. 
The so called Kuramoto-Sakaguchi model has the form 

e^u- Krsm{ei-i)- 13), 

where /3 is a phase shift of the coupling function. Low- 
dimensional dynamics of this model with Lorentzian g{uj) 
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FIG. 6. (color online) The mean field's amplitude r, the phase 
velocity Q., and the mean frequency of the ensemble f^na ver- 
sus the phase shift /3. Results from the theoretical predictions 
for the low-dimensional dynamics, Eqs. H39I - HO)) , are given 
with a solid line (black) for r, a dotted line (red) for Q, and 
a dashed line (blue) for fens- Results from numerical simula- 
tions, Eq. IT]), of the ensemble are shown with squares (black) 
for r, triangles (red) for and diamonds (blue) for fens- Pa- 
rameters: 7 = 0.5, K = a and a G [— 47r/9, 47r/9]. 

can be also obtained using the ansatz [2l| . with deriva- 
tions being similar as for the problem in Section HVAI It 

yields 

Bz K 

— + i{uj- i-f)z + —{ia^e'^ - ~z*e-'^) = 0, 

which was also obtained in [20]. From here, one imme- 
diately obtains the long-term evolution of -0, i.e. the 
frequency of the TW 

— = = c2) + if sin/3 -7tan/3. (39) 
at 

For the mean ensemble frequency in this model, Eq. (jl7p 
can be analytically solved and with further substituting 
stationary values for r and f2, it transforms to 

fens = w - si gn{Vl)^ tan P- 

This directly leads to 

fens + tan (40) 

which can also be straightforwardly obtained from 
Eq. (|23)) after applying the residue theorem for the in- 
tegral over LO- The stationary values of the macroscopic 
parameters of this system are illustrated in Fig. [51 There, 
it is also obvious that \fens\ < \^\ and that the frequen- 
cies have odd symmetry, features which are a direct con- 
sequence of Eqs. §^ and (^0]) . 

It is worth noting that by taking the derivative over /3 
from expressions (|39[ HQ)) , extreme values for both mean 
frequencies follow as 

Pniax/min = ± arCCOS 2^ / K , 

for il and 

Pmax/niin = ± arCCOS \/ 'J / K , 



for fens- Of course, since the population should be coher- 
ent, parameters have to be chosen such that K > Ke = 
27/ cos /3 holds. Thus, the values of the phase shift /3 that 
produce maximum deviation in both mean ensemble and 
mean field frequencies are directly obtained. 



V. SUMMARY 

The TW state in the KM, occurs whenever the symme- 
try in either the natural frequencies or the coupling func- 
tion itself is broken. For asymmetric coupling strengths 
on the other hand, the wave occurs only for certain cou- 
pling parameters. The results of this study indicate that 
in populations with an asymmetric, bell-shaped distribu- 
tion of the non-synchronized oscillators and equal cou- 
plings, the mean ensemble frequency is always the mean 
of the natural frequencies. In contrast, the TW of the 
synchronized oscillators - the mean field frequency - has 
a value between the mean and the mode of those. When 
the asymmetry originates from the coupling strengths or 
the coupling function, both mean frequencies have non- 
trivial values. In particular they differ from the mean or 
the modes of the frequency distribution. 

The case of the asymmetric unimodal frequency dis- 
tribution has quite a straightforward explanation. Here 
the frequency of synchronization is a result of the inter- 
play between the cluster of the locked oscillators and the 
whole ensemble. Synchronized oscillators tend for fre- 
quency locking at the mode, where the density of similar 
oscillators is the highest, whilst the influence of all oscil- 
lators would be balanced if it is at the mean of all natu- 
ral frequencies. The compromise is achieved through the 
self-consistent system for the phase balance and Eq. ([32)1 
follows directly from this interplay. This balance would 
be destroyed if the nucleus is either on the peak or on 
the mean. 

By increasing the number of locked oscillators, the 
mean field frequency shifts towards the less skewed side 
of the distribution, i.e. towards the mean, because the 
number and the influence of the drifting oscillators de- 
creases. Finally, once all the oscillators get locked, the 
frequency of the entrainment will be exactly the mean of 
the distribution. In this way, the influence of all oscilla- 
tors is equal. The same explanation can be applied for 
bimodal or multimodal cases, since again there is only 
one cluster of synchronized oscillators in the TW state. 

On the other hand, in the case of contrarians and con- 
formists, the physical explanation of the mean ensemble 
frequency, the frequency of entrainment and their dif- 
ference, is not yet very clear. Previous work went as 
far as showing that in some parameter spaces for oppo- 
site sign couplings, both peaks in the distribution of the 
phases lose their stability. Instead, they start to chase 
each other and are no longer a distance of tt apart, hence 
producing the TW. Here we have emphasized that this 
behavior simply acts to preserve the phase balance, which 
directly gives the value of the entrainment frequency. 
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while its influence on arranging the non-synchronized os- 
cillators also makes it responsible for the mean ensemble 
frequency. Despite this, we believe that even by show- 
ing, analytically and numerically, that the velocity of the 
mean phase and the frequency of the ensemble differ in 
this case, an important characteristic that seems to be 
unique for TW states only is immediately revealed. This 
explanation was missing when the model was first in- 
troduced 0, and clear distinction between macroscopic 
frequencies was not made. Moreover, "the mean phase 
velocity" obtained through the low-dimensional dynam- 
ics was defined as the mean ensemble frequency. 

In the more general and realistic situation, with asym- 
metric frequencies and distributed couplings, the macro- 
scopic frequencies can have values in very unexpected 
boundaries - even outside the limits of the natural fre- 
quencies. The reason for this complex behavior is again 
a result of the system's self-arranging, such that the in- 
fluence of the non-synchronized oscillators to the mean 
fleld amplitude vanishes and the phase balance is main- 
tained. In order for this to be achieved for multimodally 
distributed couplings, there are different clusters of syn- 
chronized oscillators for each mode. This case also shows 
possibilities for further research either in eventual hys- 
teresis dynamics and seeking the number of stable so- 
lutions, or in the clustering phenomena for distributed 
coupling strengths. 

When the asymmetry is induced through the coupling 
function itself the macroscopic frequencies always differ, 
in a similar way to the case when it stems from the cou- 
pling parameters. In that sense the mean ensemble fre- 
quency is of the same sign, but with lower absolute value 
than the mean fleld frequency. Also, their dependance on 
the phase shift is nonlinear, and the ways they respond 
to the shift do not coincide, with both of them having 
the highest values for different phase shifts. 

Taken together, these results suggest that whenever 
the population is experiencing a TW, the locking of the 
oscillators is at a frequency different to the mean of the 
instantaneous frequencies. Since in inverse problems or 
in experiments it is often only the macroscopic param- 
eters that can be obtained, a clear interpretation of the 
observed mean frequency is always needed. Moreover, 

I 



the asymmetric scenarios tend to be far more abundant 
in the real physical systems. That is to say, asymme- 
try of the natural frequencies can immediately make any 
scenario more realistic, whilst few of the examples with 
opposite coupling strengths can be traced in inhibi tory 
and excitatory neurons [233 O'" social dynamics [18| . 
As for the phase-shifted coupling function, although it 
was introduced to describe the formation of nonlinear 
waves in non-oscillatory media [ll|, it can be used to 
model various phenomena, such as Josephson junctions 
Q , mammalian intestines and heart cells . Hence, in 
the models based on measurements, like those describ- 
ing the brain dynamics Q , one should precisely define to 
which of the macroscopic frequencies, of the mean field 
or mean ensemble, the measured frequency corresponds. 
Finally, this work reveals the strong need for future re- 
search that should explain the physical link between the 
observed mean frequency of any system with cooperative 
dynamics and the two macroscopic frequency parameters 
described here for non-trivial cases. 
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Appendix A: General formula for the mean 
ensemble frequency 



We work in the reference frame rotating with 57 = t/i. 
In this frame, the locked oscillators are frozen and one 
has to consider the drifting ones only, Eq. ([T^ . This can 
be applied in the expression for the mean frequency of 
the ensemble, Eq. (ITlT) . which becomes 



fens ^ I n {uj - Krsm9)g{iu)T{K) y}'^^\ ^^''^l dwdRdO ^ h - h- (Al) 

I-ttJ J\uj\>\Kr\ 2Tr\Lj - Kr sm{e)\ 



r 



The second integral requires more calculation. Integrat- 
ing over the phases first, for positive frequencies u > \Kr\ 
ujT{K)g{uj)dKduj, (A2) it yields 

'\uj\>\Kr\ 

because of the probability normalization ^27r 



The first integral is 
h - 



27r 



2Tr\uj — KrsinB 



■de ^ 1. 



(A3) 



- {KrYKr sin 6* 
27r|w - Kr sin6l| 



(A4) 
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while for negative frequencies uj < —\Kr 

r^"" - {Kr)^Kr sin 

Jq 27r|a; - Kr sine] 

,(A5) 

It is interesting to note that both integrals are even in K. 
Following this, the second integral, I2, simply becomes 



d0 = uj + yju? - {Krf 



-\Kr\ 



-!{u})T{K)[uj + yjuj^ - {KrY]dujdK 



— 00 J —00 
00 /'OO 



g{oj)T{K)[u - v/w2 _ {^Krf]dujdK. (A6) 



00 J I Kr 



After partial cancelation of I2 with Ji because of (|Aip 
one obtains a simple formula for calculation of the mean 
ensemble frequency in the reference frame of the order 
parameter 



Jens — 



00 /"OO 



Kr \ . 

00 /"OO 



g{Lo)V{K)y/uj'^ - [KrfdujdK - 



-\Kr\ 



g{uj)T{K)y/Lo^ - [KrfdujdK 



Kr\ 



[g(u) - g{-w)]V{K).JZF~{KrfdujdK. 



(A7) 



Note that g{uj) is centered in — and hence is asymmetric, even for symmetric g{uj) which is often the case. Returning 
to the original frequencies uj, the last expression becomes 



00 poo 



\Kr\ 



[g{uj + g{^LO + VL)]T{K)yJuj'^ - {KrfdujdK. 



(A8) 
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